arXiv: 1501.0092 lv2 [math.PR] 25 Aug 2015 


PHASE TRANSITION FOR A CONTACT PROCESS WITH 
RANDOM SLOWDOWNS 

KEVIN KUOCH 


Abstract. Motivated by a model of an area-wide integrated pest management, we 
develop an interacting particle system evolving in a random environment. It is a gen¬ 
eralized contact process in which the birth rate takes two possible values, determined 
either by a dynamic or a static random environment. Our goal is to understand the 
phase diagram of both models by identifying the mechanisms that permit coexistence 
or extinction of the process. 


1. Introduction 

In this paper, we consider a generalized contact process describing the evolution on a 
lattice of a multitype population. Developed in the hfties by E. Knipling and R. Bush- 
land (see [IZI[I9]), the Sterile insect technique emerged to control the New World screw 
worm, a pest threatening livestocks in America. It is a pest control method whereby 
sterile individuals of the population, to either regulate or eradicate, are released. While 
sterile males compete with wild males, they eventually mate with (wild) females so 
that offsprings reach neither the adult phase nor sexual maturity, reducing the next 
generation. By repeated releases, we should be able to cause a variety of outcomes 
ranging from reduction to extinction. 

During the last decades, mosquito-borne diseases have become a serious worldwide 
spreading disease affecting millions of people each year. As the main vector of dengue 
fever, the mosquito-type Aedes aegypti, adapted its lifestyle from the forests to urban 
areas across the world, it is considered as one of the most fastest growing disease. At 
short notice with no available cure, population control is the only viable option to stem 
the disease. 

Female mosquitoes taking a blood meal from a virus-infected person become carriers 
of the virus. After some days, the virus spreads to the mosquito’s salivary glands 
and is therefore released into its saliva. From then on, this mosquito remains lifelong 
infected and transmits the virus to each new person it bites. As of today no cure exists, 
tackling the dengue transmission involves eradicating its vector: the mosquito. Lauded 
for several aspects and its success, the sterile insect technique is one of today main 
actors in this hght, as recent programs launched by Brazilian and European institutions 
highlight. 

To understand this phenomenon, we construct a stochastic spatial model on a lattice 
in which the evolution of the wild population is governed by a contact process whose 
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growth rate is slowed down in presence of sterile individuals, shaping a dynamic random 
environment. 

The contact process with random slowdowns (CPRS) {rit)t>o has state space {0,1, 2, 3}^ . 
Each site of Z'^ is either empty (state 0), occupied by wild individuals only (state 1), by 
sterile individuals only (state 2), or by both wild and sterile individuals (state 3). We 
only consider the type of the individuals present on each site (and not their number), it 
is biologically reasonable to assume enough females are around, not limiting matings. 
Sterile males are spontaneously released at rate r. The rate with which wild individuals 
give birth (to wild individuals) on neighbouring sites is Ai at sites in state 1 and A 2 at 
sites in state 3. We assume that 

A 2 < Ai (1.1) 

to highlight a competition between individuals occurs at sites in state 3 so that the 
fertility is decreased. Deaths of each population occur at all sites at rate 1, they 
are mutually independent. With no further restriction, this process is called symmetric 
while in a so-called asymmetric case, sterile individuals stem births on sites they occupy. 

A number of models with dynamic random environment has been investigated (see 
Broman [ 5 ], Remenik [ 21 ], Steif-Warfheimer [ 32 ], Garet-Marchand m)- On the other 
hand, identifying mechanisms yielding extinction or coexistence of species from ecology 
issues has been devoted in several papers (cf. Schinazi [31], Durrett-Neuhauser [TT] 1 
using multitype contact processes. The common feature is a type of individuals evolving 
as a contact process coexisting along with an other type evolving independently in a 
markovian way and dictating somehow its growth. 

The question we address now is for which values of r does the wild population survive 
or die out. Along with the monotonicity of the survival probability, we prove the 
existence and uniqueness of a phase transition with respect to the release rate r for 
fixed growth rates Ai and A 2 . 

Denote by Xc{d) the critical value of the basic contact process on Z'^. By comparing 
the process with basic contact processes, we explore in Propositions 12.31 and 12.41 varying 
growth parameters where no phase transition occurs i.e. when both Ai and A 2 are 
smaller or larger than Xc{d). The most interesting cases are discussed in the following 
results. 

Theorem 1.1. Suppose X 2 < Xc{d) < Ai fixed. Consider the symmetric CPRS. There 
exists a unique critical value rc G (0, cxo) such that the wild population survives ifrKVc 
and dies out if r > r^. 

Theorem 1.2. Suppose Xfid) < Ai fixed. Consider the asymmetric CPRS. There exists 
a unique critical value Sc G (0, 00 ) such that the wild population survives if r < Sc and 
dies out if r > Sc. 

These results agree with the former conclusions done by E. Knipling (1955) in a 
deterministic model. As a consequence, we discuss in Subsection 15.11 the competitive 
ability of the sterile individuals: we prove the critical value increases as the competi¬ 
tiveness of the sterilized population decreases or as the fitness of the wild population 
increases. 

To prove our results, we first construct the graphical representation of the process in 
Section [3l setting a percolation structure on which to define the process, lending itself 
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to the use of percolation that we shall do in Section |5l Afterwards, we point out general 
properties of the system in Section 01 such as necessary and sufficient conditions for the 
process to be monotone, then, only sufficient conditions in Propositions 14.21 and 14.31 to 
be in line with the construction of the process. It follows that the survival probability 
is monotone, see Corollary 14.61 The tricky part to prove these conditions lies in the 
dehnition of an order on the state space {0,1, 2, 3}^'^, since an element of {0,1, 2, 3} 
does not correspond to the number of particles but a type. 

Then, we use block constructions (see 0 ) and dynamic renormalization techniques 
(see e.g. [2]) to prove Theorems 11.11 and 11.21 in Section [5l 


So far, we were unable to get a hand on bounds for the critical rate. To this purpose, 
we hrst consider the mean-held equations in Section El featuring the densities of each 
type of individuals, we can explicitly hnd equilibria yielding numerical bounds on the 
phase transition. We derive a rigorous proof of the convergence of the empirical densities 
to these equations in [22] . 


In a very different view, we consider the process (^t,a;)t>o, where is a contact 

process on (0,1}^ with varying growth rate determined by a quenched random envi¬ 
ronment u initially hxed. In words, it describes an initial release of sterile individuals 
within the target area and looking at the induced time-evolution of the wild population. 
Using former results obtained by T.M. Liggett [23ll2l|, we obtain in Section [71 several 
survival and extinction conditions for the process. In that way, we consider two kinds 
of growth rates in Z: one where the rates depend on the edges and one where the rates 
depend on the vertices. This yields numerical bounds on the phase transition for the 
process to survive or die out. 

Random parameters for the contact process have become a matter of interest by 
showing up unexpected behaviours and interesting survival conditions of the process, 
see Andjel [1], Bramson-Durrett-Schonmann [8], Klein [18], Newman-Volchan [28] . 


The paper is organized as follows. We begin by introducing the model and results for 
a phase transition in Section [21 In this view, we provide a graphical construction of the 
process in Section El and stochastic order properties in Section 01 Section El is devoted 
to prove the phase transition on Z'^. After what, we investigate the mean held-model in 
Section El and the quenched model in Section [71 both exhibiting bounds on the critical 
value. 


2. Settings and results 

2.1. The model. On the state space O = , where F = {0,1, 2, 3} and S = Z'^, 

the contact process with random slowdowns (GPRS) is an interacting particle system 
whose conhguration at time t is rjt E O, that is, for all x G Z'^, r]t{x) G F 
represents the state of site x at time t. Two sites x and y are nearest neighbours on 
7/ if ||a: — y\\ = 1, also written x y, and ni{x,rit) stands for the number of nearest 
neighbours of x in state i, i = 1,3. We dehne now the symmetric and asymmetric 
GPRS we mentioned earlier as follows. 
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0 ^ 1 

at rate Xini[x,rj) + X 2 nfix,rj) 

1 ^ 0 

0^2 

at rate r 

2^0 

1 ^ 3 

at rate r 

3 ^ 1 



3^2 


Definition 2.1 (Asymmetric CPRS). The asymmetric CPRS {rit)t>o has transition 
rates in x for a current configuration r] : 


( 2 , 1 ) 


Definition 2.2 (Symmetric CPRS). The symmetric CPRS {rit)t>o has transition rates 
in X for a current configuration rj given by dZD, to which one adds the following 
transition: 

2 —)■ 3 at rate Xinfix, r}) + \ 2 nfix, rj). (2.2) 

Therefore, the 2’s dictate independently the growth rate of the process, and even 
inhibit births in the asymmetric case, so that the 2’s shape a dynamic random environ¬ 
ment for the I’s. 

In both cases, it rj E fl and x G denote hj rjl E fl, i E {0,1, 2, 3}, the configura¬ 
tion obtained from r] after a flip of x to state i: 

rjiu) it u ^ X 


T] —> rfi at rate c{x,rifi), where Vm G rjKu) = 


it u = X 


Let C be the infinitesimal generator of {r]t)t>o, for ciny cylinder function / on fl, 

3 

^fiv) = - fiv)) 

i=0 

with infinitesimal transition rates, common to both cases (see (M), 

Xini{x,r]) + X 2 nfix,r]) if r]{x) = 0 


(2.3) 


(2.4) 


c{x, r/, 0) = 1 if ri{x) E {1, 2}, c(a;, rj, 1) = 


c{x,ri,2) = 


r if ri{x) = 0 
1 if rilx) = 3 ’ 


1 if ri{x) = 3 


c{x, ri,3) = r if r]{x) = 1 


(2.5) 


and add the following rate in the symmetric case: 

c(x, r], 3) = Xinfix, r]) X 2 nfix, rj) if r]{x) = 2. 

Since all the rates in fl2.5p are bounded, by [26l Theorem 3.9] there exists a unique 
Markov process associated to the generator fl2.4p . Denote by {rif)t>o the process starting 
from A, i.e. such that rjo = 1 a, in other words r]o is the configuration with sites in state 
1 in A and empty otherwise. We care about the evolution of the wild population, i.e. 
individuals contained in sites in state 1 and 3. Define 

= ( 2 . 6 ) 

as the set of sites containing the wild population at time f > 0. 

Denote by PAi,A 2 ,r the distribution of (?7t°^)t>o with parameters (Ai, A 2 ,r). When Ai 
and A 2 are fixed, we simply denote this law by instead of PAi,A 2 ,r- The process {r]t)t>o 
with initial configuration 770 = l{o} is said to survive if 

> 0 , / 0 ) > 0 


(2.7) 
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and die out otherwise. 

Remark when r = 0, the process (r/t)t>o is a basic contact process, denoted by 
whose generator, dednced from (12.dh with r = 0, acts on local fnnctions of hi so that 

= E {/«3 - /«) + A.(/«i) - /K))}. 

xeZ‘^ 

It is known the basic contact process exhibits a nniqne phase transition on whose 
critical valne is denoted by Xc{d), see IHl [25l [26] for fnrther details. 

Dehne the critical value according to the parameter r by 

Tc = rc(Ai, Aa) := inf{r > 0 : P,(3f > 0, iff ^ = 0) = 1} (2.8) 

Indeed, the class {0, 2} is a trap: as soon as Ht = 0, the wild popnlation is extinct 
while sterile individnals (2) are infinitely often dropped along the time. 

2.2. Results. As annonnced in Section [H we begin by a first set of conditions for the 
process to snrvive or die ont, when Aa < Ai are both smaller or larger than Ac: 

Proposition 2.3. Suppose Aa < Ai < Xc{d). For all r > 0, both symmetric and 
asymmetric GPRS with parameters (Ai,Aa,r) die out. 

Proposition 2.4. Suppose Xc{d) < Aa < Ai. For all r > 0, the symmetric GPRS with 
parameters (Ai,Aa,r) survives. 

We now tnrn to Theorems o and II.2[ for which we shall prove: 

Theorem 2.5. Assume Aa < Ac < Ai fixed. Let {rit)t>o he the symmetric GPRS. Then, 
(i) there exists vq G (0, oo) such that if r < vq then the process {rjt)t>Q survives. 

(a) there exists ri G (0, oo) such that if r > ri then the process {r]t)t>o dies out. 

Theorem 2.6. Assume Xfid) < Ai and Aa < Ai fixed. Let {r]t)t>o the asymmetric 
GPRS. Then, 

(i) there exists Sq G (0, oo) such that if r < sq then the process {rjt)t>Q survives. 

(a) there exists Si G (0, oo) such that if r > si then the process {r]t)t>o dies out. 

Following Theorems 12.51 and 12.61 monotonicity argnments (given by Corollary 14.bh 
yield the existence of the critical valne Tc fl2.8p . We complete the stndy by investigating 
the behavionr of the critical process. 

Theorem 2.7. In both cases, the critical GPRS dies out: Pcc(Vf > 0, 7^ 0) = 0. 

We shall prove these resnlts in Section O 

3. Graphical representation 

In view of nsing percolation, we constrnct the process from a collection of independent 
Poisson processes, see [15]. Think of the diagram x M+. For each x G Tfi, consider 
the arrival times of mntnally independent families of Poisson processes: {A^ : n > 1} 
with rate r, : n > 1} and : n > 1} with rate 1 and for any y snch that 

y ^ X, {T®’^ : n > 1} with rate Ai. Let {Uf : n > 1} be independent nniform random 
variables on (0,1), independent of the Poisson processes. 
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At space-time point (x, A^), put a to indicate, if x is occupied by type-1 indi¬ 
viduals (resp. empty), that it turns into state 3 (resp. state 2). At (x, (resp. 

at put an “X” (resp.“o”) to indicate at x, that a death of type-1 (resp. of 

type-2) occurs. At times draw an arrow from x to y and two kinds of actions oc¬ 
cur according to the occupation at x: if x is occupied by type-1 individuals, the arrow 
indicates a birth in ?/ of a type-1 individual if y is empty or in state 2 ; if a; is occupied 
by type-3 individuals giving birth at rate A 2 < Ai, check at (x,T*’^) if < A 2 /A 1 to 
indicate, if success, that the arrow is effective so that a birth in ?/ of a type-1 individual 
occurs if y is empty or in state 2. In the asymmetric case, births occur only if y is 
empty. 

See Figured] for an example of the time-evolution of both processes starting from an 
identical initial conhguration, in the space-time picture Z x R_|_. 

For s < t, there exists an active path from (x, s) to {y, t) in Z*^ x M_|_ is there exists 
a sequence of times s = sq < si < ... < Sn-i < Sn = t and a sequence of corresponding 
spatial locations x = xo,xi, ...,Xn = y such that: 

i. for i = 1,..., n — 1, vertical segments {xi] x (sj, Sj+i) do not contain any X’s. 

ii. for i = 1, there is an arrow from Xi-i to Xi at times s* and if Xj_i x Sj is 
lastly preceded by a “ •” this arrow exists only if Usi~^ < A 2 /A 1 . 

and in the asymmetric case, substitute ii. by 

ii’. for i = 1, there is an arrow from Xj_i to Xi at times Si while {xi} x Si is 
not lastly preceded by a while if Xj_i x Si is lastly preceded by a this 
arrow is effective if Usf~^ < A 2 /A 1 . 

Consider the process {Af^)t>o, the set of sites at time t reached by active paths starting 
from an initial conhguration Aq = A, containing sites in state 1 in A and 0 otherwise: 

A^ = { 1 / G Z : 3x G A such that (x, 0) {y, t)} 

Then Af = with Hf dehned in fl2.6p so that A^ represents the wild population at 
time t starting from an initial conhguration A of type-1 individuals. 

This graphical representation allows us to couple CPRS starting from diherent initial 
conhgurations by imposing the evolution to obey to the same Poisson processes. Other 
kinds of couplings are possible through the analytical construction of the process as 
we will see. More generally, graphical representations allow to couple processes with 
diherent dynamics as well, we investigate this in the next section. 

4. Stochastic order and attractiveness 

The stochastic order between processes is related to the total order dehned on the 
set of values F = {0,1,2,3}. In a biological context, setting an order between types 
of individuals does not make sense, but mathematically it allows us to construct a 
monotone model and to compare diherent dynamics. To avoid confusion, denote by A 
the state 2, by B the state 0, by C the state 3 and by D the state 1. For two given 
conhgurations G fl, we have if coordinate-wise 

Tj^^\x) < ri^^\x) for all x G Z'’*, 
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time 



Figure 1. Starting from r]Q = l{o}, following the arrows, if 17° < ^ and 
U 2 < the wild popnlation occnpies at time t the set Ht = {—1,0,1} 
in the asymmetric case and the set = {—2, —1, 0,1} in the symmetric 
case. 


with respect to the following order on F\ 

A<B <C <D (4.1) 

It is important to note that this order is the only one possible with which such a 
stochastic order for the the process {rit)t>o holds. 

We now settle necessary and sufficient conditions, then only sufficient, to obtain 
properties of stochastic order which give monotonicity properties and comparisons be¬ 
tween processes. Using Borrello’s results [6], let x,y G 7A be two neighbouring sites 
and a, /5 G F to rewrite the transition rates for k G (1, 2} thanks to 

• is the birth rate of a type-1 individual in y when r]{y) = (3 and depending 
only on the value of r] in x given by a. The state in y flips from [3 io (3 + k. 

• is the jump rate of a site in state (3 to state f3 + k. 

• P~^ is the jump rate of a site in state a to state a — k for k < a. 

Theorem 4.1. |6l Theorem 2.4] For all {a, (3) G F"^, (7,5) G such that {a, (3) < 
(7,5) (coordinate-wise i.e. a < 7 and (3 < S), hi > 0, ji > 0, an interacting particle 
systems (ylt)t>o with transition rates (F°’^, F“^) is stochastically larger than an 
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interacting particle system {Bt)t>o with transition rates Pa */ 



where and = P-\ 


One has for the asymmetric GPRS, the following rates. 



(4.3) 


to which, one adds the following rates if we consider the symmetric GPRS, 



(4.4) 


(4.5) 


(4.6) 


Anticipating forthcoming compntations, we expand fl4.2p when {At)t>o and {Bt)t>o 
correspond to two symmetric GPRS with respective rates (Ai, A 2 , r) denoted with upper 
indexes (2) and (1). Note that by fl4.3p - f|4.6p . the following inequalities can be used as 
well for asymmetric GPRS or basic contact processes. Taking ji, hi G {0,1} (one can 
check they are the only non trivial possible values), we have respectively from fl4.2p fii- 
(ii) applied to our model: 


l{ji = 0}l{fc = 2}1{S -P= l}(l{5 = G, d = B}{R°jfPh{a = D} + R°S%^^h{a = C}) 

+ 1{6 = B,P = A}{R°jffh{a = D} + R°S%Pl{a = G})) 

+ l{ji = 0}l{fc = 2}1{S - d = 0}(l{<5 = d = = D} + R°Sm^h{a = G}) 

+ 1{5 = I3 = A}{R°^'^fh{a = D} + R°sT^l{a = G})) 

+ l{ji = 0}l{k = 1}1{S - d = 0} (l{<5 = d = + 1{<5 = d = A}pP^^) 

+ l{ji = l}l{fc = 2}l{<5-d = 0} 

(l{S = d = = D} + R°S%^^h{a = C}) 

+ 1{S = (3 = A}{R°p,^fh{a = D} + R°S%^^h{a = G})) 

< l{ji = 0}1{Z = 2}(l{^ = = D} + R°S%^^h{^ = G}) 

+ 1{6 = A}(i?^"f R{7 = D} + <"/>1{7 = G})) 

+ l{ji = 0}1{; = 1} (l{(5 = + 1{S = A}pP^^) (4.7) 

+ l{ji = 1}1{1 = 2}(l{S = B}{R°^P^h{y = D} + R°S%^^h{y = G}) 
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and 

l{hi = 0}l{k = l}(l{a = + l{a = 

+ l{hi = 0}l{k = 2}(l{a = D}P-'^’^^^ + l{a = 

+ l{hi = l}l{k = 2}(l{a = + l{a = 

> l{hi = 0}1{Z = 2}1{7- a = l}(l{7 = D,a = + 1(7 = C,a = 

+ l{hi^ 0}1{Z = 2}1{7 - a = 0} (l{7 = a = D}P~^’^^^ + 1(7 = a = ClPJ 
+ l{/ii = 0}l{Z = l}l{7-a = 0}(l{7 = a = L»}P"^’^^^ + l{7 = a = C'}P“^’^^^) ^ 

+ l{hi = 1}1{; = 2}1{7 -a = 0}(l{7 = a = PjP”^’^^^ + l{a = 7 = 

We are now ready to state and prove the following resnlts. 

Proposition 4.2. The symmetric CPRS is monotone, in the sense that one can con¬ 
struct on a same probability space two symmetric CPRS cif^d with 

respective parameters and (Af\A^^^ such that 

^ Vt^'^ C: a-s. for all t > 0 (4.9) 

if and only if all parameters satisfy 

1 . A« < aS'\ 5. aS'^ < A®, 5. r(^) > r( 2 ) 7 , yW < 

2 . \f> <\f\ I A<‘> < Af, «■ aS” < 1. «. r<‘l > 1. 

Remark conditions 1. and 2. are given by assnmption (ll.lh . 

Proof Let and be two symmetric CPRS with respective parameters 

(A^^\ A 2 ^\r^^^) and (A[^\ By Theorem 14.11 necessary and snfficient conditions 

on the rates for to be stochastically larger than are given with {a, ft) < 

( 7 ,( 5 ), by (14.71) - (14.8p . All different possible scenarios provide the following necessary 
conditions: 

(I) ji e {0, l},5 = (3 e {A, B} in (gTD give 

(i) (T = C, 7 = D: A 2 ^^ < Af\ This is a conseqnence of conditions 1. and 3. 
or 2. and 4. 

(ii) a = j = C\ stated by condition 4. 

(iii) a = 'y = D: A^^^ < A® stated by condition 3. 

(II) 71 = 0. d = ^. (5 = 1 + d = C in m) give 

(i) a = D: A[^^ < 1 stated by condition 6 . 

(ii) a = C: A^^^ < 1 stated by condition 7. 

(III) hi = 0, 7 = q; e {B, D} in fl4.8p give r^^'> > stated by condition 5. 

(IV) hi = 0,a = B,'y = l + a = C in (14. 8 p give r^^'> > 1 stated by condition 8 . 

while in other scenarios, one retrieves rednndantly the above conditions or tantological 
ineqnalities snch as “1 > 0”. Finally, we obtained the conditions stated 1. to 8 . □ 
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Following Proposition 14.21 we now constrnct a coupled process on x 

such that According to the given order fl4.1|) on F, as 

ni{x,7]^Q^) +n^{x,r]^Q^) < + n^{x,r]^Q^). 

In what follows, we construct the coupling through generators. The three following 
tables depict the inhnitesimal transitions of the coupled process. 



Table 1. Coupled transitions of two symmetric CPRS, via a basic cou¬ 


pling, that preserve the partial order. 


transition 

rate 

(0,3)^ < 

r ( 1 , 1 ) 
( 0 , 1 ) 
( 2 , 2 ) 

1 (2,3) 

(1 - Aj^^)ni(x,7?P)) + (1 - A^^^)n3(x,7?P)) 

1 

7 .( 1 ) _ 1 


Table 2. Coupled transitions of two symmetric CPRS, not via a basic 
coupling, that preserve the partial order. 
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Table 3. Coupled transitions of two symmetric CPRS, via a basic cou¬ 


pling, that do not preserve the partial order. Nevertheless, starting from 
an initial conhguration that does, the coupled dynamics does not reach 
these states. 


To verify all the rates above are well dehned, one decomposes and 

713 ( 0 ;,(i = 1 , 2 ), as follows 

ni{x,r]^^'^) = \{yr^x: = 1 }| 

+ \{y^x-. r]^^\y) = = 3}| + |{7/ ~ x : ri^^\y) = e {0,2}}|, 

n 3 {x,r]^^^) = \{yr^x: y^^\y) = v^^\y) = 3}| + |{7/ ~ x : y^^\y) = e {0,2}}|, 

ni(x, = \{yr^x: y^^\y) = V^^\y) = 1 }| 

7i3(x,77(^)) = \{yr^x: y^‘^\y) = r]^^\y) = 3}| + \{y ^ x : y^‘^\y) = l,v^^\y) = 3}|, 

in which case, we decompose the rate 

A®ni(x, 77 *-^^) — A^^^ni(x, 77 *-^^) -|- A®7i3(x, — X^2 ^n3{x, 

= (Af^ - a5^^)|{7/ ~ X : y^^\y) = y^^\y) = 1}| + (A^^ - a5^^)|{7/ ~ x : y^‘^\y) = l,y^^\y) = 3}| 

+ ~ ® = 3}| +>^?\{y ~ x : y^‘^\y) = l,y^^Hy) e {0,2}}| 

+ ~ X : y^‘^\y) = 3,y^^\y) € {0,2}}| (4.10) 

which is non-negative under conditions 1. to 4. coming from (I) and (III) in inequalities 

(HZD-dH). 
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Rates of Table |2] are non-negative thanks to conditions 6. to 8., given by inequalities 
(II)-(i)(ii) with (3 = B,6 = C. Condition 5. is used by Tables [U and |3] that correspond 
to a basic coupling while Table [2] uses a different coupling. Table |3] listing transitions 
of the coupled process starting from conhgurations that do not preserve the dehned 
partial order, nevertheless, starting from an initial conhguration that does, dynamics 
of the coupling given by Tables [U and [2] do not reach states of Table [3] so that marginals 
are retrieved. 

For a coupled process starting from an initial conhguration such that 

^ since transitions of the two hrst Tables preserve the order on F, the mar- 
kovian coupling we just constructed is increasing: 

pl'Jo = 1 for alH > 0 (4-11) 

where stands for the distribution of {ri^\ starting from {r]^\ 

We now investigate sufficient conditions for attractiveness, under which we build a 
coupled process through a basic coupling as follows. 

Proposition 4.3. The symmetric process {r]t)t>o is monotone, in the sense that one 
can construct on a same probability space two symmetric processes {ril^'^)t>o o,nd 

with respective parameters and{Xf\ X^\r^‘^^) satisfying G {0,1}^'^, 

such that 

for all t>0 a.s. (4.12) 

if the parameters satisfy 


1 . 5. AS')<Af\ 

9 \(2) / \(2) , x(l) / x(2) 

/Cf. A2 ^ , 4 . A2 ^ A2 , 


3. ^ 


Remark 4.4. In view of the proof of Proposition \4 ■ 3[ one can actually relax the admis- 
sible initial conditions: it is enough to assume that and satisfy rjj^^ < rjl^^ and 
for all X G {ril^\x), rij^\x)) 7 ^ (0,3). In particular one could start from rjl^^ = rj^\ 

Proof. Given initial conditions, possible states for the coupled process keep laying in 
Table [Hof Proposition 14.21 and the coupled process does not reach any state of Tables [2] 
and 121 One can therefore omit conditions 4. to 6. of the previous Proposition 14.21 and 
transition rates from the couple (0, 3) can be dehned through a basic coupling even if 
it does not preserve the order: 


transition 

rate 

(0,3)—^ < 

r (1,3) 
(0,1) 
(0,2) 

1 (2,3) 

A(^^ni(a:,?7p)) + X^ 

1 

1 

^(1) 


Table 4 


in which case. Table IH substitutes Table El 

If (hi^^)t>o and {'riP)t>o diher by at most one parameter, one deduces 


□ 
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Corollary 4.5. Suppose € {0,1}^'^. Then for the processes {ril^^)t>o (I'^d 

with parameters and (Af\A^^^ respectively, one has 

(i) Attractiveness: z/(A[^\ A^^^ = (Af\A^^^ , then ==» 17 ^^^ < 

a.s., for allt > 0. 

(a) Increase w.r.t. Ai; if{pt^^)t>o o,nd{pf'*)t>o have respective parameters {X^i \ X 2 ,r) 
and (A[^\A 2 ,r) such that X 2 < A^^^ < A^^, then p^^ < pj^^ pj;^'' < pf'^ a.s., 
for all t > 0. 

(Hi) Increase w.r.t. X 2 : if{Pt^^)t>o o^nd {pf'*)t>o have respective parameters {Xi, X^\r) 
and {Xi,Xf\r) such that A^^^ < A^^^ < Ai, then p^^ < pj^^ pj:^^ < pf'^ a.s., 
for all t > 0. 

(iv) Decrease w.r.t. r: if{pi^^)t>o o.nd {pl'^^)t>o have respective parameters {Xi, X 2 ,r^^'>) 
and (Ai, A 2 , such that r^^'> > r^'^l with X 2 < Xi, then Pq^'^ < pj^'' pf^^ < pf"* 
a.s., for all t > 0. 

A highlighted consequence related to Corollary (I4.5p -(iv) is the monotonicity of the 
survival probability with respect to the release rate r for fixed Ai, A 2 : 

Corollary 4.6. Suppose X 2 and Ai fixed. If {pt)t>o has initial configuration po = l{o}, 
the mapping 

r ^ P,(Vf > 0 , Ht^ 0 ) 

is a non-increasing function. 

Proof. Let {pl^'’)t>o and ip?^)t>o be two processes such that p^\p^^ G {0,1}^'*, with 
respective parameters (Ai, A 2 ,r‘^^^) and (Ai, A 2 ,r‘^^^) such that r^) < r^‘^\ then by Corol¬ 
lary |4]5l for alH > 0, 

ffP) (- ffW ^ fjm ^ 

□ 

Analogously to Propositions 14.21 and l4.3l we can obtain conditions for the asymmetric 
CPRS to be attractive. For the rest of the section, we refer the reader to [21] for more 
detailed proofs. 

Proposition 4.7. The asymmetric CPRS {pt)t>o is monotone, in the sense that con¬ 
ditions 

i. A«<A«, S, AS‘'<Af>, 5. r-m^rW, 

9 \P) ^ \P) ! \<‘l ^ \P) 

A2 ^ A]^ ; 4. A2 ^ A2 ; 

are sufficient to construct on a same probability space two asymmetric CPRS (? 7 i^^)i>o 
and {pt^^)t>o with respective parameters A^^^ and (Af\A^^^ and with 
initial condition p^\p^^ G {0, such that 

for all t>0. (4-13) 

Proof. As in the proof of Proposition 14.31 one applies P, Theorem 2.4] with ji,hi G 
{0,1} to two asymmetric processes {p^'^)t>o and {pf^)t>o with respective parameters 
(Ai^\ A 2 ^\ rd)) and (Af\A^^\ r^)). Using fl4.7p - (14.81) with rates given by (14.31) . we de¬ 
duce the following necessary conditions from the first inequality: 
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(I) ji e { 0 , 1 }, 5 = /S = 5 in (021) give 

(i) a = 'y = C, [] = B, 6 = C: stated by condition 4. 

(ii) a = C, 7 = iA: A^^^ < A[^^ stated by conditions 1. and 3. 

(iii) a = 7 = ZA; A^^^ < A^^ stated by condition 3. 

(II) ji = 0, /3 = S, 5 = 1 + /9 = C in (14.Zp give 

(i) a = D: aS^^ < 1. 

(ii) a = C: A^^^ < 1. 

The relation (14.Sh staying unchanged, one has 

(III) hi = 0, 7 = a G {B, D} in fl4.8p give stated by condition 5. 

(IV) hi = 0,0; = i ?,7 = l + o = C'in (14. 8 h give > 1. 

To sum up, the obtained necessary conditions are 


1. a" < a", 3. A« < Aft 5. r<‘> > r'^) 

2. Af> < Af, 4. Af> < Aft 6, Af> < 1. 


7. A^^^ < 1, 

8 . > 1 . 


As for Proposition I4.2[ these conditions allow us to construct an increasing markovian 
coupling and are thus necessary and sufficient. But as in Proposition 14.31 given our 
initial conhgurations, state (0, 3) is not possible for the coupled process. One can thus 
dispense conditions 6 to 8. and conditions 1. to 5. are sufficient to construct a basic 
coupling. □ 

We now compare asymmetric, symmetric and basic contact processes. 

Proposition 4.8. Let {r]t)t>o be an asymmetric GPRS and {xt)t>o be a symmetric 
GPRS, both with parameters (Ai,A 2 ,r) and r]o,xo G {0,1}^'^ such that A 2 < Ai, then 
for all t > 0, 

ho < Xo ^ ht < Xt a.s. for allt>0 

Proof. Consider a symmetric process with rates fl4.3p - fl4.4p and an asymmetric process 
with rates (I4.3p . Similarly to the proof of Proposition 14.21 we apply inequalities (14.7p - 
(14. 8 p to these rates so that they give the following necessary conditions: 

(I) ji e {0, 1}, 5 = jd = B,a = C, 7 = ZA in (14.7p give: A 2 < Ai 
(II) hi = 0, a = ZA, 7 = 1 + a = C in (14. 8 p give r > 1 
As previously, condition r > 1 is necessary to construct an increasing markovian coupled 
process in a general framework, but if one restricts the initial conditions to satisfy 
ho ^ Xo and rjo, xo G {0,1}^ , condition A 2 < Ai is sufficient and the coupled process 
can be constructed through a basic coupling. □ 

Proposition 4.9. Let {^t)t>o be a basic contact process on {0, with growth rate 
Ai and let (xt)t>o be a symmetric GPRS with parameters (Ai, A 2 ,r) such that A 2 < Ai. 
Then, 

Xo<io^Xt< it a.s. for all t > 0 

Proof. Consider a symmetric CPRS with rates fl4.3p - fl4.4p and a supercritical contact 
process {it)t>o viewed as a symmetric CPRS with parameters (A[^\ \ with A^^ = 

Ai, A 2 ^^ = 0, = 0, that is with rates fl4.5p . Applying fl4.7p - fl4.8p to rates fl4.3p - fl4.4p 

and fl4.5p . since values A and C do not exist for the process ft) we have the following 
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necessary condition: ji G {0, 1 }, /9 = 5 = i?, a = C, 7 = D in (14.7p give A 2 < Ai. 
While relation (I4.8p does not give further condition. Condition A 2 < Ai is sufficient and 
allows us to construct a coupled process through basic coupling. □ 

Proposition 4.10. Let {rit)t>o be a symmetric CPRS with parameters (Ai,A 2 ,r) and 
(6)t>o o basic contact process on {2,3}^^* with rate A 2 such that A 2 < Ai. Then, 

Co < I/O ^ 6 < ht a.s., for all f > 0. 

Proof. For the process (Ct)t >05 values B and D are not reached so that the necessary 
and sufficient conditions on the parameters given by fl4.7p - f|4.8p exhibit the condition: 
ji G {0,1}, /3 = (5 = y4, a = C, 7 = H in fl4.7p give A 2 < Ai. Inequality fl4.8p gives 
no condition on the rates and condition A 2 < Ai is sufficient to construct the coupled 
process {^t,Vt)t>o via a basic coupling. □ 

5. Phase transition on Z"* 

Here, we take advantage of all the stochastic order relations established in Section 0] 
to derive a phase transition for the CPRS, in both symmetric and asymmetric cases. 
According to fl2.7p . assume the process to have initial conhguration rjQ = l{o} and note 

T]t = We first explore cases when A 2 < Ai are both smaller or larger than Ac. 


Proof of Proposition 12.31 Let (^t)t>o be a basic contact process with growth rate Ai 
and let (7t)t>o be a symmetric CPRS with parameters (Ai, A 2 ,r) such that ijo < ■Co- By 
Proposition I4.91 is stochastically larger than {rjt)i>Q. Since Ai < Ac, is 

subcritical, thus, the symmetric CPRS dies out. 

The extinction of the asymmetric CPRS is a consequence of the extinction of the 
symmetric process and Proposition I4.8I □ 


Proof of Proposition \2.4\ Let be a basic contact process with growth rate A 2 on 

{2, 3}^ and let (i/Joo be a symmetric CPRS with parameters (Ai, A 2 ,i’). By Proposi¬ 
tion ITTUl {^t)t>o is stochastically lower than {r]t)t>o- Since A 2 > Ac, the process 
is supercritical and therefore, the symmetric CPRS survives. □ 

Before proving Theorem I2.5I to I2.7I , the following Subsection I5.1I deals with conse¬ 
quences of Theorems I2.5I and I2.6I along with monotonicity results of Section 01 

5.1. Behaviour of the critical value with varying growth rates. Here we prove 
monotonicity between growth rates and the release rate in the sense that, the more 
competitive the species is (i.e. the higher the parameter A 2 is) or the fitter the species 
is (i.e. the higher the parameter Ai is), the higher the release rate is (i.e. the higher 
the critical value Vc is). 

Proposition 5.1. For j = 1,2, the function Xj \— > i'c(Aj) is non-increasing. 

Proof. We consider j = 2, the case j = 1 is similar. Let {rjt)t>Q and (7t)t>o be two CPRS 
with respective parameters (Ai, A2,r) and (Ai, A2,'r) such that A2 < A2. Existence of Tc 
and for those processes is given by Theorems 12.51 and 12.61 By contradiction, suppose 
Tc > r(, 

Vt "L Vt by Theorem 12.71 and Corollary 14.51 fivL 

P,(Vf > 0, Hi ^ 0) < P,/(Vf > 0, Hi ^ 0) = 0 


Let r be fixed be such that > r > rl. By Corollary 14.51 fiiiL if 70 = Vo then 
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But since r < Vc, the process (?7t)t>o survives: Pr(Vt >0, Ht ^ 0) > 0. This contradicts 
Vt < Vt hence Tc < r^. □ 

5.2. Subcritical case. The following proof relies on a comparison with an oriented 
percolation process on the grid C = {{x,n) G : x + n is even,n > 0}. We follow 
arguments used in [20] for a spatial epidemic model. 

Proof of Theorem A2.3 (i). To simplify notations, choose d = 1 but the proof remains 
the same for any d > 2. Introduce the space-time regions, 

B = (—4L,4L) X [0,T], Bm,n = {2mLei,nT) + B 
I = [-T, L], Im = 2mLei + I 

for positive integers L,T to be chosen later, where (ei,...,ed) denotes the canonical 
basis of 


Consider the process restricted to the region Bm,n, that is, constructed from 

the graphical representation where only arrival times of the Poisson processes occurring 
within Bm,n are taken into account. Therefore, a birth on a site x G Bm,n from some 
site y only occurs if y G Bm,n- By Proposition 14.31 and Remark [4.41 one has 


m,n 

Vt 


< Vt 


(5.1) 


for alH > 0 if = ho L 

Let k = \ \fL\ and define C = [—k, k]. One declares {m,n) G £ to be wet if for any 
conhguration at time nT such that there is a translate of C full with I’s in with 
containing only O’s and I’s, the process restricted to Bm,n is such that at time (n-f 1)T 
there are a translate of C in Im-i and a translate of C in Im+i, both full of I’s, with 
Im-i and Im+i containing only O’s and I’s. 


We show the probability of a site (m, n) G £ to be wet can be made arbitrarily close 
to 1 for L and T chosen sufficiently large. By translation invariance, it is enough to 
deal with the case {m,n) = (0,0). 


Suppose / contains only O’s and I’s as well as the translate of C full of I’s and set 
r = 0 in 13, that is, no more type-2 individuals arrive in the box B after time 0. If type-2 
individuals are present on (—4L, —L) U (L,4£) x {0}, the probability of the event E 
“they all die by time T/2” is at least 

(l-exp(-T/2))'^ 

which is larger than 1 — e for T and L large enough. On £, the process restricted 
to the box B is now from time T/2 a. supercritical contact process (^™’"')t>r /2 with 
distribution G •). Define the hitting time of the trap state 0 of the contact 

process starting from [—£, P\ and restricted to [—£, £] x [0, T/2]. To ensure there are still 
enough I’s for Mountford has proven 

P{T{i) < exp(£)) < exp{—i) for £ large enough (5.2) 

Partition C into M = \ \/k\ boxes, each of them being a translate of [0, Mj. From each 
of these M boxes, say box j < M, run a supercritical contact process (Ct )t>o which 
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coincides with the restriction of to this box. Then for the union of these j boxes 
the probability there is at least M I’s within C by time T/2 is after fl5.2p . with T such 
that exp(M) < T/2, at least 

P(r(M) > T/2)^ > P(r(M) > exp(M))^ > (1 - exp(-M))^^ (5.3) 

which can be made larger than 1 — e, for M, i.e. L, large enough. 

R. Durrett and R. Schinazi [12] have shown that for a contact process (^t)t>o, for any 
A C Z, except for a set with exponentially small probability, either Ef = 0, or 
on a linearly time growing set [—at, at]: there exists a > 0 such that for all R C Z, 
there exist constants C, 7 such that for a; G A + at, 

P(S,'' 0, 7(x) ^ < Cexp(- 7 t) (5.4) 

We applied this result with A G C, the set of I’s in the box. We just proved that 
|x4| > \M\. Moreover by attractiveness (i.e. Proposition 14.2p . one can choose k, and so 
L, large enough so that this supercritical contact process starting from at least M 
I’s survives at time T/2 with probability close to 1, hence 7^ 0 and fl5.4p is valid. 
In this situation, taking T/2 = 9L/{2a) with L large enough, the process (^P’"')t>o 
starting from at least M I’s in [—T, T] at time T/2 will be coupled with a process on 
[—3T, 3T] with probability at least 1 — e at time T. Hence, since 3T > aT/2 > 2L, by 
time T the contact process started inside [—T, T] has not reached the boundary of 
[—4T,4T] with probability close to 1. Then, the process and the contact process 
are the same with probability 1 — e in [—4T, 4T] ; this way, the coupling of (^P’"')t>o 
with (Cf)t>o works so far with probability 1 — e if T is large enough. 

Since the distribution of (Cf)t>o is stochastically larger than the upper invariant 
measure P of the contact process, on the survival event, P loads a positive density p of 
I’s. Since P is ergodic, 

1 

Z 1{4(L = 1}=P P-a.e. 

x=—3L 

Then if L is large enough, there are at least k I’s in any interval of length 2T + 1 with 
P-probability at least 1 — e. Since we obtained that (^r’")t>o is coupled to (^f)t>o by 
time T with probability at least 1 — 2e, for L large enough, there are at least k I’s in 
[—3T, —T] at time T with probability at least 1 — 2e and similarly, at least k I’s at time 
t in [L, 3T] with probability at least 1 — 2e as well for Consequently, 

P((0, 0) wet) > 1 — 4e, if r = 0. (5.5) 

Since H is a hnite space-time region, for hxed L, T, one can pick tq > 0 small enough 
so that if r < ro 

Pr((0, 0) wet) > 1 — 6e. 

Now construct a percolation process by dehning Gm,n = {(lu, n) wet}. Note Gm,n de¬ 
pends only on the process constructed in Bm,n, and for (a, b) G C, events Gm,n and Ga,b 
are independent if (m, n) and (a, b) are not neighbours. The events {Gm,n, (nr, n) G £} 
are thus 1-dependent. By the comparison theorem CQl Theorem 4.3], the process 
(7™’"')t>o restricted to regions Bm,n is stochastically larger than a 1-dependent percola¬ 
tion process with probability 1 — e. 
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One can choose e small enongh so that percolation occnrs in the 1-dependent perco¬ 
lation process with density 1 — e, see m Theorem 4.3]. □ 

5.3. Supercritical case. In the following, one compares onr particle system with a 
percolation process on x Z_|_ nsing arguments exposed as in [30] . 

Proof of Theorem \2.5\ (ii). Assume d = 2, the proof can similarly be extended to higher 
dimensions. For all {k,m,n) G Z^ x Z_|_. Introduce the following space-time regions, 
for positive L and T to be chosen later. 

A = [—2L, 2L]^ X [0,2T] Ak,m,n = A + {kL, mL, nT) 

B = [-L, Lf X [T, 2T] Bk,n^,n = B+ {kL, mL, nT) 

C ^bottom {_}^side ^k,m,n ^ T {kL,TnL,nT) 

where Cbottom = {{'m, n,t) G .4, : t = 0} and 
Cside = {{^,n,t) G A : \m\ = 2L or |?7,| = 2L} 

Consider a restriction of the process {r]t)t>o to Ak,m,n, that is, the process 
constructed from its graphical representation within Ak,m,n- 

One declares a site {k, m, n) G Z^ x Z+ to be wet if the process {rj^ ’^’^)t>o contains 
no wild individual in Bk,m,n starting from any configuration in Ck,m,n- Therefore it will 
be the same for ?7t I . . Sites that are not wet are called dry. 

For any e > 0, we show that for some chosen L and T any site {k,m,n) is wet with 
probability close to 1 when r is large enough. By translation invariance, it is enough to 
consider {k,m,n) = (0,0,0). Set r = cxd in 4,. Then, the process contains 

only sites in state 2 or 3: sites in state 0 or 1 flip instantaneously into state 2 and 3 
respectively. That is, {rit''^’'^)t>o is in fact a contact process {^’'^’^)t>o on {2, 

The contact process {f,t)t>o on {2, 3}^^ with growth rate A2 < Ac is subcritical. 

If there is some wild individual in B, it should have come from a succession of births 
started somewhere in C. Starting from a site in Cside, a path to B should last at least 
L sites ; according to C. Bezuidenhout and G. Grimmett [5] there exists such a path 
with probability at most Gexp(—7L), for some positive constants G,7. Hence, 

P(3(x,t) G Cside X [0,2T] : {x,t) B) < 4(2T x (4L-|- l))Gexp(—7L) 

Similarly, starting from the base Cbottom, there exists a path lasting at least T sites with 
probability 

F{3{x,t) G Cbottom ■ {x,t) B) < (4L l)^Gexp(-7T) 

Gonsequently if r = 00, 

P((0, 0, 0) wet )> 1 - 4(2L x (4L + l))Ge-'>'^ - (4L + > 1 - e/2, 

for L and T large enough. 

Since 4. is a finite space-time region, one can pick r large enough so that with 
probability at least 1 — e/2, an exponential clock with parameters r rings before any 
other so that there are no type-1 individuals in A with probability close to 1: 

Fr{k, m, n) wet) > 1 — e 


for r large enough. 
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To construct a percolation process on 1? x one puts an oriented arrow from 
(fc, m, n) to (x, y,z)iin < z and if Ak,m,rP^-Ax,y,z 7 ^ 0- The event Gk,m,n = {(^, in, n) wet} 
depends only on the graphical construction of the process within Ak,m,n, furthermore, 
for all {k,m,n) G x there is a hnite number of sites {x,y,z) G x Z+ such 
that Ak,m,n n Ax,y,z 7^ 0- The percolation process is dependent but of hnite range. The 
existence of a path of wild individuals for the particle system corresponds to a path of 
dry sites for the percolation and we proved that dry paths are unlikely. This implies 
that for any site, there will be no more wild individual after a hnite random time, see 
[31 Section 8 p.l6]. By translation invariance, this corresponds to the extinction of the 
process. □ 


This proves the existence of a phase transition for the symmetric CPRS. The proof 
of Theorem 12.51 fil only uses that contact process with growth rate Ai is supercritical, 
hence the existence of sq in Theorem l2.6f fil as well. By Proposition 14.81 the asymmetric 
CPRS dies out if the symmetric one does so that existence of si in Theorem l2.6R iil 
follows from Theorem l2.5R iih Though, one remarks conditions of Theorem 12.61 are 
milder: one can actually show the existence of Si in a neater way, assume A 2 > Ac: 
without transition ”2 —)• 3” in the asymmetric case and choosing r = 00 , notice for a 
subcritical contact process on {2,3}^^, there are no paths of wild individuals created 
by the 3’s from the boundary Ck,m,n up to extinction, but this occurs exponentially fast, 
see [5]. 


5.4. Critical case. Note the arguments developed by [3] rely on elementary properties 
of the contact process making them robust. They are also valid for the CPRS since it 
satishes the following properties. 

(A) contact process-like dynamics: one retrieves the growth rate Ai or A 2 of a basic 
contact process, even if it is determined randomly. 

(B) attractiveness, cf. Section HI 

(C) correlation inequalities: such as FKG inequality. 

Note that the use of (C) is possible because the proof works in hnite space-time regions 
and since survival events are monotone functions of Poisson processes (from the graph¬ 
ical representation, see Section [3]), correlation inequalities for such Poisson processes 
apply, see [T6l Lemma 2.1] for instance. We therefore omit the proof of Theorem 12.71 
but rather present a brief outline and refer the reader to |IlEaE2] for further details. 


Outline of the proof of Theorem fS. 7| The hrst step consists to observe that if the 
process survives in an arbitrary large box, then, that it reaches its boundaries densely. 
We shall estimate these densities at each side of a space-time region. 

This way, one can repeat this step by running the process in an other adjacent box 
starting from the boundary of the previous one and so on, conditionally on the fact 
that the starting conhguration is dense enough. This is the second step. In connection 
with the proofs of Theorems 12.51 and 12.61 where we looked after having translations of 
occupied hnite intervals at a given hxed time, here we look after having translations of 
the densities in some space-time slab. 

Now, knowing that at each stage, one can construct overall a path of adjacent boxes 
wherein the process survives and reaches the boundaries densely, it remains to compare 
the process with an oriented percolation process to extend the arguments to inhnite 
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space and time. As before, compare a space-time box to a vertex in the even lattice of an 
oriented percolation so that one declare a vertex to be wet if some good event associated 
to the box is a success. Conclude thanks to well-settled results about percolation theory. 

To conclude the section, Theorem ll.il (resp. Theorem 11.21) follows from Theorem 12.51 
(resp. Theorem 12.61) and by monotonicity arguments given by Corollary 14.61 

6. Mean-field model 

The mean-held model associated to the CPRS is a deterministic and non-spatial 
model where all individuals are mixed up, featuring the densities of each type of indi¬ 
viduals. Such models give rise to differential systems and are interesting to compare 
with stochastic behaviours. We investigate here the equilibria, in both asymmetric and 
symmetric models, since existence of such equilibria yield the existence of a critical 
value. We exhibit conditions on r to deduce bounds on the critical value. 

Subsequently, let Uj be the density of type-i individuals for i = 1,2,3. Overall, 
one has ui + U 2 + = 1 — uq. Furthermore, in connection with the dehnition of 

wild and sterile individuals, we consider as well ui, resp. ^2, the density of the wild 
individuals (type-1 and type-3 individuals), resp. the sterile individuals (type-2 and 
type-3 individuals), and the density of empty sites uq = uq. Relations between the 
M-system and the n-system are described by 

{ Ui = l-Vo-V2 

U2 = 1-Vq-Vi . ( 6 . 1 ) 

U:i = VQ + Vi+V2-l 

Since we consider densities, both systems satisfy 

WiG [0,1] fori = 0,1, 2, 3, G [0,1] for i = 0,1, 2. (6.2) 

Following the most suitable case, we investigate either the n-system or both systems. 

6.1. Asymmetric CPRS. Assuming total mixing, the mean-field model associated to 
the asymmetric CPRS is given by: 


v'q = -2d{^\2 - Ai)no + ^2^1 (A2 - Ai)ni + Xi - X2^vo - {r + 2 )vo - vi- V2X-2 

v'l = 2d(^{X2 — Ai)'Uo -l- A2ni -I- (A 2 — Ai)n2 -I- Ai — X2^vo — vi 
V 2 = r(l - V 2 ) — V 2 


This system gives rise to two equilibria: 

= (0,0, 


(6,3) 


r -|- 1 


(r -h 2)" 


. 2 -|-t’ r -\- 2 

vo,vi^V2 VddAir-I-2(iA2r ’ 2(r-I-1) 2(4dAi-|-2(iA2r) ’ r-|-l/ 

Note that the hrst equilibrium gives {ui,U 2 ,U 3 ) = ^0,^,0^, which puts a positive 
density on the sterile individuals and none on the others, which corresponds to the 
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extinction of the process. By checking conditions fl6.2p on the second, we get that: the 
density vi is non-negative as soon as 

4(iAi -I- 2 d\ 2 r > (r -|- l)(r -|- 2). 

which gives the following condition 

2d\2 - 3 + ^(2dA2-3)2-8(l-2dAi) 


r > 


( 6 , 4 ) 


This indicates a lower bound for the phase transition. 

6.2. Symmetric CPRS. The mean-held equations associated to the symmetric CPRS 


are : 


u{ = 2 d{\iUi + \ 2 U^)uo + U 3 - {r + 1)mi 
u'2 = ruo + U3- U2- 2d{XiUi + X 2 U 3 )u 2 
Mg = rui + 2d{XiUi + X2U^)u2 — 2uz 

As previously, this system admits one trivial equilibrium: 


(6.5) 




+ 

retrieving once again a situation related to the extinction of the process, by a positive 
density of sterile individuals and none of the wild ones. We derive the non-trivial 
equilibrium thanks to the corresponding M-system: 

v'q = -2d[{X2 - Ai)mo + A2M1 + (A2 - Ai)mi + Xi - X2^vo - {r + 2 )vo - V1V2 + 2 

v[ = 2d(^{X2 — Ai)Mo + -^ 2^1 + (-^2 ~ •^l)P 2 + Ai — A2^)(1 — Mi) — Ml 

V 2 = r(l - M 2 ) - M 2 

( 6 . 6 ) 

Let us determine the non-trivial equilibrium. Last line of fl6.6p gives already M2 = 

T 

-. Using relations of fl6.ip in fl6.6p . according to m( = 0, an equilibrium (mo,Mi,M 2) 

r + I 

satishes in particular 

Ml = 2(i(AiMiA2M3)(1 - Ml) (6.7) 

checking mi cannot be equal to 1, one then has 

Ml 


Ml 


— 2 d{XiUi -|- A2M3) 


( 6 . 8 ) 


and 


Ml 7^ 1. (6.9) 

On the other hand, from the M-system fl6.5p with relations fl6.1l) and using condition 

dnn). 


m'i = 0 + (2 r)vo + Ml 


- Ml r + 1 

, , , r -|- 2 MnMi 

M 2 = 0 (r -F 2 )mo + Vi -—r 


= 0 


Mg = 0 (1 - Mi)^ + (1 - Mi)( 


r -|- 1 1 

1 


Ml 


r -I- 1 


(r -F 1 )mo) - Mo = 0 
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By solving the last line with respect to (1 — ni), one has 


1 — ni = (r + l)no or 1 — ni = 


r + 1 


But since 1 — ui G [0,1], necessarily 1 — ui = (r + l)no- Using fl6.8p . vq solves 
2(i(Ai + A2'r)(r + l)'yg — (2(iAi + 2dX2T + r + l)no + 1 = 0. 

This implies 


no = 


or Vq = 


2dXi + 2dX2r 


r + 1 

(1) if Vq = then by relations flO.ip . 

Ml = 0, Ms = 1, M2 = 1 + M2, 


which is a contradiction. 
(2) if Vo = 777^— nrirz:’ fhen 


2dXi + 2dX2r ‘ 


Ml = 


r + 1 


V2 = 


2d{Xi + X2r) ’ r + 1’ 

Verifying this M-system to be a set of densities by condition fl6.2p . one case highlights 
a condition on r: Mi < 1 if and only if r(l — 2dX2) < 2dAi — 1. In the case where 
-^2 < l/(2(i), then one has the condition 

2dAi - 1 


r < 

- 1 - 2dX2 

Consequently, a non-trivial equilibrium of fIb.Op is given by 

( \ ‘ ^ M + 1 

[Vo,Vi,V2) = 


( 6 . 10 ) 


( 6 . 11 ) 


^2dX\ + 2dX2r^ 2dXi + 2dX2r^ r + 1/ 

To put in a nutshell, this survey of equilibria associated to both mean-held models, 
in symmetric and asymmetric case, gave us the bounds fl6.4p and fl6.10p for the phase 
transition. 


7. Quenched random environment on 

We now derive bounds on the phase transition by using T.M. Liggett’s results [251124] 
for the one-dimensional contact process in quenched environment, that we adapt to our 
different dynamics. Let p G (0,1), dehne a random environment uj G {0,1}^ where each 
site a; G Z is either free (0) with probability 1 — p or slowed-down (1) with probability 
p, independently of any other site. 

The contact process in random environment (CPRE) {^t)t>o that we consider has 
state space {0,1}^ with quenched environment cm. Let Ai and A2 be growth parameters 
such that 

A 2 < Ac < Ai, (7.1) 

where Ac is the critical rate of the one-dimensional contact process. Subsequently, 

p = r/(r + l) (7.2) 

for r G (0, cx)), is the probability a site is slowed down. Deaths occur at rate 1. If 
{u:{x)^x G Z} is chosen deterministic (hence inhomogeneous), we will refer to it as the 


















PHASE TRANSITION FOR A CONTACT PROCESS WITH RANDOM SLOWDOWNS 


23 


contact process in inhomogeneous environment (CPIE). Both CPRE and CPIE are still 
monotone according to Section 01 

Denote by the distribntion of the CPRE with parameters (Ai,A 2 ,r) and en¬ 

vironment u. For hxed parameters Ai and A 2 , simplify by P“. For any A C Z, dehne 
Ef := {x ^ = 1}, where denotes the process at time t started from the 

initial conhgnration ,^0 = 1a- If R = {0}, simplify by Et = For almost-every 

realization of u, the CPRE is said to survive if 

p;f(Vf >0, St ^ 0 ) > 0 

and die out otherwise. 

Consider snbseqnently two kinds of random environment: one depending on the 
vertices and one depending on the edges of the graph. The following resnlts are adapted 
from [23l Theorems 1 and 3] and [231 Theorem 1.4]. 

7.1. Random growth on vertices. If Ay(k) is the growth rate from site /c G Z (See 
Fignre |2]): a birth at site k occnrs at rate A^(A: — 1) if fc — 1 is occnpied pins at rate 
At,(A; -|- 1) if fc -|- 1 is occnpied, where 

Ay(k) = Ai(l — u(k)) + A2cu(k) (7.3) 


A,(A:-2) A,(A:-1) A,(k) A,(k +1) 



Theorem 7.1. Assume u is a fixed environment. The CPIE dies out if for all n G Z, 

k n 

sn Av{j + 1) < 00 and ER A^(j - 1) < 00 . (7.4) 

k>nj=n k<nj=k 

If the family {u{k), G Z} is i.i.d. then the family {Av{k), /c G Z} is i.i.d as well and 

Corollary 7.2. The CPRE dies out i/log A.u(0) < 0. That is, if A 2 < I and 
r > - log A 2 / log Ai. 


Theorem 7.3. The CPRE survives if 


j>0 


1 Ay(k) + Ay{k — 1 ) -|- 1 \ 

A„{k)A„{k-l) ) 


< 00 . 


No independence in the above prodnct disables ns to obtain explicit conditions for 
snrvival of the process. Nevertheless, by dehning the randomness on the edges rather 
than on the vertices, meaning that the growth rates emanating from a site k respectively 
to fc -|-1 and to fc — 1 are randomly chosen for each A: G Z, we are able to explicit bonnds 
on r with respect to Ai and A 2 . 
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7.2. Random growth on oriented edges. Let {(pe(fc), ^e{k)), fc G Z} be an ergodic, 
stationary and i.i.d. seqnence. Here, given a site /c G Z, a birth from k to k + 1 occnrs 
at rate Xe{k+ 1) and independently of a birth from k to k — 1 occnring at rate pe{k — 1). 
See Fignre |3l 


Xe{k-1) X^{k) Xeik + 1) Xe{k + 2) 



Figure 3. Random environment on oriented edges 


Snppose both rates are i.i.d, defined for all A; G Z by 

Xe{k + 1) = Pe{k - 1) = Ai(l - Uj{k)) + X2Uj{k). (7.5) 


Theorem 7.4. The CPRE dies out if 


¥XfXe{k) < 1 and 1 — E); 


Xe{k) 


< Et 


Xe{k) 


l-EfX,ik) 


Ai 


that is, if X 2 < I and r> 

1 — X 2 

Theorem 7.5. The CPRE survives if for all j > 0, 


e:; 


Xe{j + 1) 


e; 


, Ae(0) + Pe(0) + 1 A 
Ae(0)pe(0) 


< 1 


that is, if X 2 < I + \/2 < Ai and r < X 2 (Ai — \/2 — l) / (Ai (A2 — \/2 — l)) 


7.3. Nnmerical bonnds on the phase transition. Back to the basic contact pro¬ 
cess, assnme r = 0 then for all a; G Z, uj{x) = 0 a.s. and Xf,{x) = Pe{x) = Ai. We thns 
recover the one-dimensional basic contact process with growth rate Ai. In this case, 
onr estimates lead to the following bonnd for Ac. 


Corollary 7.6. Eor the one-dimensional basic contact process, 

Ac < 1 + \/2. 

Proof. According to Theorem 17.51 if r = 0 the process snrvives if Ai > 1 + \/2. □ 


This bonnd is qnite rongh bnt we derived it simply. Conseqnently, one first dednces 
a bonnd on the critical valne Ac on Z: Ac < 1 -l- \/2 ~ 2.41. 


From previons resnlts, let us derive numerical bounds on the segment where the phase 
transition occurs. Choosing Ai and A2 given by Theorem 17.41 gives us lower bounds on 
the phase transition and satisfying condition from Theorem 17.51 gives us upper bounds. 


Ai 

A 2 

phase transition 

1000 

0.8 

[0.49,4995) 

100 

0.8 

[0.48,495) 

10 

0.8 

[0.36,45] 

2 

0.8 

(0,5] 


Ai 

A 2 

phase transition 

1000 

1.4 

[ 1 . 37 , 00 ) 

100 

1.4 

[1.34, 00 ) 

10 

1.4 

[1.04, 00 ) 

2 

1.4 

M+ 
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Remark that the necessary condition A2 < 1 disables us to conclude to an upper 
bound for values of A2. In a similar way, Theorem 17.51 imposes Ai to be larger than 
1 + \/2, disabling us to hnd an explicit lower bound in such cases. 

7.4. Proofs. The contact process in random environment {^t)t>o on {0, 1}^ introduced 
by T.M. Liggett [2S1121] has transitions at x G Z given by 


0 —>■ 1 at rate p{x)^{x + 1) + A(x)^(x — 1) 
1 —0 at rate (5(x) 


(7.6) 


where the family p(x), A(x)), x G Z} stands for the random environment, which 

is an ergodic stationary process. To apply |23l Theorems 1 and 3] and [H] Theorem 
1.4] to our case, the link to fl7.3p and fl7.5p is given by: for all k E Z, 

X{k + 1) = p{k - 1) = Xv{k) = Pv{k) 

X{k) = Xe{k) and p{k) = Pe{k) 

Proof of Corollary \7.^ By the ergodic theorem. 



lim i V log A^(j) = E]:' log^(O), 


where Ef log.j,(0) 


rlogA2 + logAi 

-, so that by theorem 

r + 1 


, so that by Theorem 17.11 extinction occurs if 


r > -logAi/logA2. 


□ 


Proof of Theorem \7.4\ Since KfXe{k) = Ai(l — p), we have both conditions 


A2 < 1 and r > (Ai — 1)/(1 — A2). 


(7.7) 


Meanwhile, the second assertion of Theorem 17.41 holds if the polynomial 



is positive which holds if 
(Ai + A2 — 


2)(Ai + A2) + (Ai — A2 )'\/(Ai + A 2 — 2)^ + 4Ai A2 


4Ai(l — A2) 


but this condition is satished as soon as (17.7p are. 


□ 


Proof of Theorem\7 . 5\ Condition of Theorem 17.51 is satished if the following polynomial 




Since we assumed A2 < Ac, where Ac < 1 + \/2 by Corollary 17.61 only conditions of 
Theorem 17.51 remain. □ 
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